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Abstract 

Taking the same trial wave function, the ground-state energy of Frohhch po- 
laron is investigated by variational method and coherent-state expansion (CSE) 
one, respectively. Within the accuracy to a^(the electron-phonon coupling con- 
stant), both methods can give the same analytic expressions of polaron ground- 
state energy as the function of electron-phonon coupling constant. We find that 
the CSE method can much simplify the calculation and shows more advantages in 
higher order approximations. 

PACS numbers: 71.38.-k, 63.20.Kr, 31.15.-p 



Electron-phonon (e-p) interaction is an important elementary excitation in con- 
densed matter physics. Although the Green Function method is quite mature to deal 
with such problems, it still can't provide ideal results when the e-p interaction is rel- 
atively strong, such that, the variational method were often involved in consideration. 
Polaron problem is particular notable in ionic or semiconductor materials. Various ap- 
proaches had been applied to the calculation of polaron ground-state energy during 
the past 40 years. Recently, based on an effective Hamiltonian of Frohlich polaron under 
LLP transformation, a new concise method ^ was proposed by Wang et al. by ex- 
panding the wave function of phonon's coherent-state in powers of phonon's correlation. 
Putting the wave function into Schodinger equation, H| >= E\ > and comparing the 
coefficients of same power of phonon's creation operator in two sides of the equation, 
a set of coupled equations can be obtained, and rather good results of ground-state 
energy and system wave function can be calculated by iteration method, especially in 
a wide range of coupling strength materials. Differently from variational method and 
other methods, coherent-state expansion (CSE) is systematic, concise and effective ex- 
hibited in two aspects: the first is that the establishment of the wave function will be 
not changed in different problem, the second is that we can expand the wave function 
to higher order phonon's correlation. In fact, the great success of CSE method have 
been shown in the study of mobile polaronl^, exciton-phonon system^l and dissipative 
two-level system^^etc. 

The systematic CSE study with higher order phonon's correlation had been per- 
formed in ID polaron ^1 system and other systems ^HjiCiliC^l- example, on 
the physical study of two-site Holstein model^l, CSE calculations are relatively simple 
compared to other approaches. Theoretically, CSE trial wave function can be expanded 
to arbitrary orders, it is just an exact example to demonstrate the "convergence" of CSE 
method. In order to illustrate the reliability and applicability of CSE method and pro- 
vide a new choice to other physicists, we adopt the same wave function as the article jS] 
and use variational method to study the intermediate coupling polaron ground-state en- 
ergy. In 1968 Larsen had studied the polaron ground state energy and effective mass 
by using variational method with a similar variational wave function^ . He studied the 
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static polaron problem and the wave function he took is a special case of CSE method. 
In addition, the second-order explicit expression of polaron ground-state energy as a 
function of electron-phonon coupling constant had not been given in previous papers [HI, 
[S] and 0. In present paper we will discuss and compare the calculation results under 
these two methods. 

By adopting the units of 2m = uJo = V= h = l, the Frohlich polaron Hamiltonian 

reads 

H = P2 + ^ a+flq + J2 Vq(aqe'''" + a+e-'"^-'), (1) 
q q 

where 

and Oq are creation and annihilation operators of the longitudinal optical(LO) phonon 
with wave vector q, respectively, a is a dimensionless electron-phonon coupling constant, 
N is the space dimensionality. After proceeding of unitary transformation, the 

polaron Hamiltonian becomes 

H = ^(1 - 2Q ■ q + q^)ala^ + ^ ' q ^qaq'^qV + ^^^^q + "i^- 
q q,q' q 

where Q is total momentum of a polaron, Q = P + ^ qa^ Oq, and it is a conserved 

q 

quantity. In order to compare with CSE method, we assume the variational wave function 
is the same as Refs. [6] [7] 

I >= I >o + YKM'(^tK'\ >o 
q.q' 

where 



>r 



Jjexp[aqa^]|0 > . (5) 



Based on the eqOl ftq and 6q^q/ are supposed to be the real variational parameter. The 
ground-state energy of the system is 
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E = -^r:-- (6) 



The aim of present paper is to prove the rehabihty of CSE method from another 
point of view, therefore we only expand the formula El into a square power of e-p coupling 
constant and the terms with higher order a are neglected during the derivation. One 
finally has 
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E = ^ Aqtt^ + 2Vqaq + J] q ■ q'c 
q q,q' 

+ 2 ^(&q,q'(v4q + v4q/ + 2q ■ q')«qaq/ + VqOq/ + V^q/Ctq) 

q,q' 

+ 2^(62,^,(Aq + Aq, + 2q-q')) (7) 
q,q' 

where Aq = 1 — 2Q ■ q + q^. By variation energy ((Zj) with respect to Oq and feq,q', we 
get 



(Aq + Aq, + 2q ■ qOttqCtq' + VqC^q. + Vq^ttq + 2(6q,q,(Aq + A^, + 2q • q') = (9) 

It is easy to get the iterating equations as follows 

aqftq/ Vqftq/ + Vq/ ftq 

2 2(Aq + Aq, + 2q-qO ^'^^ 



q - --^ + 4- y2 I '^q'(^q + A' + q • q')"q + Sttqaq/l/q/ + aq,Vq + Vq 
^q ^q , I 



V^qCtq/ + V^q/Ctq 

ttn = 7^ + — > <^ «;/(An + An/ + q • q + /Ittnan'l/n' + + Vn'-j — T 

/Iq + /Iq/ + 2q ■ q 

(11) 

Taking the initial value of Oq = — Vq/Aq, substituting it into (fTU)) and (fTTj) to 
proceed the iteration, and then put it into (|7j), thus the exact expression of ground-state 
energy as the function of e-p coupling constant can be derived as following 

q ^A2Aq,(Aq + ylq,+2q-q')- ^ ^ 



The second energy term of N-dimensional polaron reads 

Taking N=3 and Q = 0, the system ground-state energy can be given by 
E = -a- 0.0159196a2 (14) 
For CSE methodic 1^, the ground-state energy and coefficient 6q^q/ are given by 
E = Q^ + J2v^a^, (15) 



- -Aq + Aq, + 2q.q' 

In terms of Eqs. (fT3|) - (fT7|) . the same expressions of ground-state energy (fT^ - (fT^ 
can be obtained when it is taken within the accuracy of a^. 

From above calculations we see, the variational method gives the same analytic 
expression of polaron ground-state energy with CSE method in the accuracy of a^. This 
energy formula had been proved correct by comparing it to the exact numerical value 
in jH]. Because the variational method is a very mature technique, our calculation may 
help us to understand the CSE method better in another point of view. 

As an example, if the trial wave function is taken in a third-order approximations, 
the variational derivation will show more complex than CSE method and yield a tedious 
mutually coupled equations. But for CSE method, it can easily get the third-order(see 
(fT ^ .(fTB j) .(fT7 |) ) ground state energyl^, which can be compared nicely with the sixth 
perturbation method results by Smondyrev'^ and SelyuinEl. It can be concluded that 
the CSE is a good approximate method. 

We should emphasize that, although the same analytic expressions of ground state 
energy can be obtained in the same accuracy of by two methods, they are intrinsically 
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different in the calculation of higher order approximations. Comparatively, CSE method 
has a more precise form and a much simple calculation. We beheve that CSE method 
will show a potential applications in other analogous systems. 
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